Griffiths-like phase in the three-dimensional ± J Edwards- Anderson spin-glass model 
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Using information of the ground state topology and the damage spreading technique we show 
a Griffiths-like phase is present in the three-dimensional ±J Edwards-Anderson spin glass model. 
With spin-flipping dynamics damage spreads for temperatures larger than T g , the glass transition 
temperature. With spin-orienting dynamics and for temperatures in T g < T < Td, damage spreads 
over a finite region of the system, composed of finite clusters of ferromagnetic character. Td is the 
spin-orienting damage critical temperature, which is of the same order as the critical temperature 
of the ferromagnetic Ising model, T c . The present results allow us to identify the clusters which 
originate the Griffiths-like phase in the present model. 

PACS numbers: 75.10.Nr, 75.40.Gb, 75.40.Mg 



The nature of the low-temperature phase of spin 
glasses has been an elusive issue over many years de- 
spite intensive studies p|. It is now well established 
that there exists a spin glass phase below a character- 
istic temperature T g for Ising spin glasses, and that far 
above T g the system is in a paramagnetic state. The 
critical temperature T g has been numerically calculated 
using diverse methods, most of them adopting the idea 
of a global growing characteristic length in the spin glass 
phase (see Ref. [2[ and references therein). However, 
it is not clear what is the origin of the 'order' behind 
the growing length. What is the underlying structure 
of the low-temperature spin glass phase also remains an 
open subject. A different approach attempted to iden- 
tify local structures presenting typical features of ordered 
states. These ideas have been further extended in the last 
decade through a systematic study of spatial and dynam- 
ical heterogeneities in spin glasses @, 0, [f| . In particular, 
the origin of strong dynamical heterogeneities is unveiled 
by ground state (GS) spatial heterogeneities in the ±J 
Edwards- Anderson (EA) spin glass model @, 0|- 

The slow relaxation properties observed above T„ have 
also attracted much attention. Randeria et al. |8j dis- 
cussed the high temperature relaxation of spin glasses. 
They found that, for models with a bounded bond dis- 
tribution, for temperatures between T g and T c , the crit- 
ical temperature of the orderer system, the time correla- 
tion function has a lower bound with a relaxation slower 
than exponential but faster than a power-law. This be- 
havior was intuitively associated with a "Griffiths-like" 
phase. Let us recall that the Griffits phase was origi- 
nally proposed for the site diluted Ising ferromagnetic 
model, where a fraction 1 — p of the sites are removed 
In the phase diagram of this system the ferromagnetic 
phase is present for temperatures smaller that T c (p) up 
to the critical percolation value p c . It was argued that 
for temperatures in the range T c (p) < T < T c ( p=1 ), sin- 



gularities in the zero field susceptibility and anomalous 
relaxation are observed due to the presence of finite size 
ferromagnetic clusters. Randeria et al. Q did not iden- 
tify unfrustratcd cluster corresponding to the Griffiths 
phase. Also, further numerical evidence does not agree 
with the predicted functional form analytically obtained 
considering Griffiths singularities It has also been 

argued that these singularities are too weak to be ob- 
served in classical systems Extensions to quantum 



models, where the sing ularities are stronger, have been 
intensively studied [12(. 

Further information on the dynamical properties above 
T g has been obtained using the damage-spreading tech- 
nique [l3j |. In this case, one measures the configura- 
tional distance between two replicas of a system, which 
are initially set to different initial conditions, and evolve 
under the same thermal noise. Damage-spreading with 
spin-orienting dynamics in the three-dimensional (3D) 
±J EA spin glass model was initially studied by Der- 
rida and Weisbuch They found three temperature 

regimes: (i) a high-temperature regime T > T\ ~ 4.1 
where the damage is zero, (ii) an intermediate regime 
T-2 ~ 1.8 < T < T\ where the damage is nonzero 
and independent of the initial condition, and (iii) a low- 
temperature regime T < T2 where the damage depends 
on the initial condition. These authors stated that the de- 
termination of T\ and Ti could be consistent with T\ = T c 
and T2 = T g . Note that this was a surprising result since 
there is no reason a priori why an imprint of ferromag- 
netic order in the dynamical properties of this spin glass 
model should be expected. 

Throughout this work, by using the damage-spreading 
technique with both spin-flipping and spin-orienting dy- 
namical rules, and its correlation with the GS topological 
properties, we studied the dynamical properties of the 3D 
± J EA model. This enables us to show that there is an 
ordered structure which presents ferromagnetic-like clus- 



tcrs capable of explaining both the anomalous relaxation 
and the ferromagnetic signature on the intermediate tem- 
perature range, T g < T < T c . Thus, these results give a 
natural framework to justify the presence of a Griffiths- 
like phase in the 3D ±J EA model. We also stress here 
that our aim is not to list the dynamical properties of 
this model with the damage-spreading technique, which 
have been extensively analyzed in the past. Instead, we 
seek for the constrained structure of the GS which un- 
veils the Griffiths-like phase in the present model, and 
thus permits us to understand all the previous puzzling 
results within a single framework. 

We consider, as a starting point, the topological prop- 
erties of the GS configurations of the ± J EA model [15j, 
E3, 0, Ell. A systematic study has shown that 
it is possible to find a backbone, characterized by sol- 
idary spins, which maintain their relative orientation 
over all the set of multiple-degenerate GS configura- 
tions [l7, 2, 2(J. In 3D, the subset of solidary spins com- 
prises a fraction of approximately 76% of the total spins. 
The cluster size distribution presents the same character- 
istics than in the diluted problem, i.e. a power-law de- 
cay form [2(J. Also, the largest cluster of the backbone 
percolates through the system. Besides, the backbone 
has a very small frustration (of the order of 10% of the 
bonds), which suggests that ferromagnetic-like order can 
be sustained at low temperatures [20|. Therefore, these 
characteristics naturally point to the backbone as the re- 
sponsible for the existence of the Griffiths-like phase in 
spin glasses in the range T g < T < T c . For tempera- 
tures below T g a ferromagnetic-like order grows within 
the largest cluster of the backbone 2l| . Another impor- 
tant characteristic is the presence of non-solidary spins, 
defined as those outside the backbone. These spins re- 
main in a paramagnetic state for all temperatures. 

In the following we will show how damage-spreading 
technique can be used to unveil the limits of the Griffiths- 
like phase. It is well known that for the pure ferromag- 
netic model the results of this technique depend on the 
Monte Carlo dynamics [bj ] . For spin-orienting dynamics 
the initial damage is healed above a damage temperature 
Xd = 4.3 [H, [23||, and it spreads to a finite fraction of 
the sample below Td- The opposite scenario is observed 
when using a spin-flipping dynamics. The damage tem- 
perature Td obtained with each dynamics is of the same 
order, and slightly lower than the critical temperature 



T c = 4.515 [24[. Figure [T] shows a schematic view of a 



± J EA system split up into solidary (green region) and 
non-solidary (violet region) spins. It is also shown how 
we expect it would respond to the location of an initial 
damage at a finite temperature T < T g = 1.120] when 
a spin-flipping dynamics is considered. Figures [TJa-c) 
qualitatively represent snapshots of the evolution of the 
system when the initial damage (red dot) is in the set 
of solidary spins, with the dashed red regions indicating 
damage spreading. Since the damaged region sustains 
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FIG. 1: (Color online) Schematic representation of spin- 
flipping damage spreading (dashed-red) below T c and con- 
sidering the separation on solidary (green region) and non- 
solidary (violet region) spins. Each set of snapshots corre- 
sponds to the evolution of the system when the initial damage 
(red dot) is within: (a-c) the set of solidary spins, (d-f) an 
island of non-solidary spins, (g-i) the set of solidary spins but 
close to an island of non-solidary spins. 



a ferromagnetic-like order, the damage should rapidly 
heals. Instead, if the initial damage is located within an 
island of the non-percolating set of non-solidary spins, as 
represented in Figs. [TJ d-f), we expect damage to spread 
over the entire island reaching a finite fraction of the sys- 
tem. An interesting particular case is when the initial 
damage is within the set of solidary spins, but close to 
an island of non-solidary spins. In this case a fluctuation 
can allow the damage to reach the border and spread 
over all the island, as represented in Fig. [ljg-i). When 
T > T g , since T is above the ordering temperature of the 
largest component of the backbone, a large fraction of 
the system will be in a paramagnetic state and the small 
initial damage will propagate to the whole system. 

The situation is different for spin-orienting dynamics. 
Above T c all the system is paramagnetic and the dam- 
age does not propagate. For T g < T < T c the largest 
cluster of the backbone is above its ordering tempera- 
ture and thus the contribution of ferromagnetic-like or- 
der comes only from the finite clusters of solidary spins, 
that are characterized by a power-law size distribution. 
Thus, these are the clusters over which damage propa- 
gates. Therefore, the damage-spreading technique over 
the finite temperature range T g < T < T c is expected to 
uncover the presence of the unfrustratcd clusters which 
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have the structure to support a Griffiths-like phase. In 
the following, through exhaustive numerical simulations 
of the damage spreading process, we shall show quanti- 
tative results that support this picture. 

The 3D ± J EA spin glass model is described through 
the Hamiltonian H = J2{i j) JijCiPj, where the sum runs 
over the nearest neighbors of a cubic lattice with linear 
size L and N — L 3 spins, <Xj = ±1 are spin variables, and 
the random bonds Jy = ± J are chosen from a symmetric 
bimodal distribution. The damage spreading protocol is 
as follows: for a given disorder realization we build two 
replicas of the system, a and /3, one resulting from equi- 
libration at a given temperature and the other from the 
same configuration but with a single damaged spin, i.e. 
with a randomly selected spin that is reversed. Equilib- 
rium configurations were obtained with a parallel temper- 
ing Monte Carlo algorithm 25[ . The two replicas of the 
system evolve under the same thermal noise, i.e. using 
the same sequence of random numbers, and the Ham- 
ming distance between each sample is measured. Thus, 
the total damage at a given time t is defined as 



(i) 



where the angular brackets and the overline correspond 
to averages over thermal histories and disorder realiza- 
tions, respectively. In this way the damage can take val- 
ues in the range < D(t) < 1. In the following we 
will present data for cubic geometry with L = 8, where 
reliable information of GS configurations can be ob- 
tained [2(| ■ The disorder average is taken over 1000 sam- 
ples, for which we have determined the backbone struc- 
ture through extensive numerical simulations. Whenever 
possible we will compare with a larger system size in or- 
der to show that our data does not suffer from strong 
finite size effects. 

We begin our analysis considering the spin-orienting 
heat-bath rule. At high temperatures the total damage 
for the spin glass model vanishes at long times, while 
at low temperatures it arrives at a stationary regime. 
Figure [^a) shows the time evolution of the total damage 
at T = 2.5 < T d for L = 8 (closed circles) and L = 20 
(open circles). The initial damage is D(0) = l/N, which 
corresponds to a flip of a random site. Note that for 
both system sizes the damage tends to almost the same 
stationary value, showing that there are no important 
finite size effects. 

Now, for a given time evolution of the whole system, 
we also perform a constrained measure of the damage re- 
stricting the sum in Eq. ^ to solidary spins, denoted 
as D s , or to non-solidary spins, denoted as D ns . No- 
tice that the normalization factor and the initial damage 
will be thus given by the average number of solidary or 
non-solidary spins, 1/N S and 1/N ns , respectively. Fur- 
thermore, as one can exactly locate the initial damage, it 
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FIG. 2: (Color online) Evolution of the damage with (a) spin- 
orienting dynamics at T = 2.5 and (b) spin-flipping dynamics 
at T = 0.6, comparing the total damage D with D^ s . The 
total damage for L = 20 is also shown. The inset shows data 
in linear scale for the stationary regime. 



is possible to have different constrained measures such as: 
(i) the damage in solidary spins when the initial damage 
is in any random spin, which we shall denote as D^, (ii) 
the damage in non-solidary spins when the initial damage 
is in any random spin, D^ s , (iii) the damage in solidary 
spins when the initial damage is within the set of solidary 
spins, , (iv) the damage in non-solidary spins when the 
initial damage is within an island of non-solidary spins, 
D^s , etc. In Fig. [2ja) we show the constrained measure 
D^s ■ The inset shows the other constrained measures as 
indicated. It can be clearly observed that the only curve 
which is appreciably different in the stationary regime is 
■D„ s s , which is smaller than the other curves. Instead, 
the damage can spread over the finite clusters of solidary 
spins below Td, resulting in D > D„ s s . 

We also consider the spin-flipping Metropolis rule. Fig- 
ure EJb) shows the total damage with this dynamics at 
T = 0.6 < T g , for L = 8 and L = 20, and also its 
comparison with -D„ s s . After a transient regime all the 
curves tend to stationary values. The inset shows the 
other constrained measures. With spin-flipping dynam- 
ics damage cannot spread over the backbone of solidary 
spins since they are far below their critical temperature, 
as depicted in Fig. [1] However, damage spreading can 
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FIG. 3: (Color online) (a) Temperature dependence of the 
damage for spin-flipping and spin-orienting dynamics. The 
total damage for L = 20 is also shown. The inset shows 
a comparison between L = 8 and L — 20 for spin-flipping 
dynamics in an enlarged scale, (b) Normalized distances be- 
tween D and -Dj^f (see the text). 



be observed for non-solidary spins. One can thus observe 
that D* s > Df, i.e. no matter where the initial damage 
is, damage spreads mostly over non-solidary spins. The 
difference is even larger between D„ S S and D^. 

The temperature dependence of the stationary damage 
is presented in Fig.[3Ja). We present the behavior of the 
total damage for L = 8 and L = 20 and also D^ s s for 
L = 8. Note that the temperature scales used for spin- 
flipping or spin-orienting dynamics are different. Both for 
spin-orienting and spin-flipping dynamics, the measure 
of the damage spreading restricted to non-solidary spins 
when the initial damage is in the subset of non-solidary 
spins, -Dn S s , gives a distinctive feature of the non-trivial 
separation originated in the GS topology of the system. 
This can be quantified by a normalized distance between 
DT and D defined as 



d s f = 



D NS -D 

D NS ' 



(L r) 



D 
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where d s t and d so correspond to spin-flipping and spin- 
orienting dynamics, respectively. The temperature de- 
pendence of d so and d s i is presented in Fig. EJb) , together 
with the relevant temperature scales, T g , Td and T c . 

We show here that spin-orienting and spin-flipping 
damage-spreading processes in spin glasses do not have 



the same damage spreading temperature as observed in 
ferromagnets. Indeed, in Fig. [3jb) one can observe that 
d s { is larger than zero below a temperature of the order 
of T g . This result indicates the ferromagnetic-like char- 
acter of the largest cluster of the backbone below T g . On 
the other hand, d so signals an important difference be- 
tween solidary and non-solidary spins close to the Ising 
damage spreading temperature Td. This difference is di- 
rectly related to the ferromagnetic character of the finite 
sized clusters of solidary spins. Indeed, this gives a reli- 
able picture allowing to understand the puzzling question 
on the origin of both Griffiths-like anomalies and ferro- 
magnetic signatures for T g < T < T c @, [H E3, [H, [27 1 . 
In fact, the GS information, used to identify the clus- 
ters which originate the Griffiths- like phase in the present 
classical problem, gives an important insight to rational- 
ize the presence of rare clusters related to the origin of 
the Griffiths phase in the quantum transverse-field Ising 
spin glass. 
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